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It is shown that natural systems are undeniable subject to random fluctuations, arising from 
either environmental variability or internal effects. In this paper, we present a spatial version of 
phytoplankton-zooplankton model that includes some important factors such as external periodic 
forces, noise, and diffusion processes. The spatially extended phytoplankton-zooplankton system 
is from the original study by Scheffer [M Scheffer, Fish and nutrients interplay determines algal 
biomass: a minimal model, Oikos 62 (1991) 271-282]. Our results show that the spatially extended 
system exhibit a resonant patterns and frequency- locking phenomena. The system also shows that 
the noise and the external periodic forces play a constructive role in the Scheffer 's model: first, the 
noise can enhance the oscillation of phytoplankton species' density and format a large clusters in the 
space when the noise intensity is within certain interval. Second, the external periodic forces can 
induce 4:1 and 1:1 frequency- locking and spatially homogeneous oscillation phenomena to appear. 
Finally, the resonant patterns are observed in the system when the spatial noises and external 
periodic forces are both turned on. Moreover, we found that the 4:1 frequency- locking transform 
into 1:1 frequency-locking when the noise intensity increased. In addition to elucidating our results 
outside the domain of Turing instability, we provide further analysis of Turing linear stability with 
the help of the numerical calculation by using the Maple software. Significantly, oscillations are 
enhanced in the system when the noise term presents. These results indicate that the oceanic 
plankton bloom may partly due to interplay between the stochastic factors and external forces 
instead of deterministic factors. These results also may help us to understand the effects arising 
from undeniable subject to random fluctuations in oceanic plankton bloom. 



PACS numbers: 87.23.Cc, 05.40.-a, 82.40.Ck 
Contents 

I. Introduction [l| 

II. Model i 

III. Results i 

A. Dynamics under the presence of noise only 

B. Dynamics under the presence of the external 
forces only [sl 

C. Dynamics under both noise and external forced 

IV. Conclusion and Discussion 
Acknowledgments d 

A. Stability analysis with the help of Maple 

References [TO 



* Electronic address: liuqx315@sina.corn| 
^Electronic address: bai-lian.li@ucr.edu 

■1- Corresponding author; Electronic address: jinzhn@263.netj 



I. INTRODUCTION 

Many mechanisms of the spatiotemporal variability of 
natural plankton populations are not known yet. Pro- 
nounced physical patterns, such as thermoclines, up- 
welling, fronts, and eddies, often set the frame for the 
biological process. While many parameters and variables 
can be learned from standard oceanographic measure- 
ments of temperature, salinity, nutrients, and biomass 
concentrations of phytoplankton and zooplankton, now 
new technologies such as remote sensing are being de- 
veloped for observing the ecosystem. Measurements 
of the underwater light field are made with state-of- 
the-art instruments and used to calculate concentra- 
tions of phytoplankton biomass, such as chlorophyll, 
as well as other forms of organic matter. More re- 
cently, satellite remote sensing and detailed numerical 
simulations identify some spatial patterns including fil- 
aments, irregular patches, sharp gradients, and other 
complex structures involving a wide range of spatial 
scales and time scales by utilizing the species' concen- 
tration Figure [1] shows two pictures deriv- 
ing from the field observation by the satellites remote 
sensing, where color gives us very useful information 
and ideas about the changes of chlorophyll concentra- 
tions in the two-dimensional spaces. The more pictures 
about the spatial patterns of the phytoplankton can be 
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FIG. 1: (Color online) This satellite images are the field ob- 
servation of the phytoplankton blooms in the Malvinas cur- 
rent region. The enhanced natural color images show ac- 
tual differences in water color while the pseudocolor image 
shows chlorophyll concentration. The images are taken form 
[http://oceancolor.gsfc.nasa.gov, with permission from Janet 
W. Campbeh.^ 



obtained from the web http: / / oceancolor.gsfc.nasa.gov" 
such as the stripelike, spotted-like, clockwise-rotating (or 
counterclockwise-rotating) spiral waves, spatial chaotic 
patterns, etc. Especially, In Ref. [1] from an observing 
system for monitoring the ecosystem is approached by 
posing five questions: How is the ecosystem changing? 
What are the forcing factors causing it to change? How 
does the ecosystem respond to natural and human forces? 
Can we predict future changes? And what are the con- 
sequences for stakeholders in our region? 

Learning from historical records and from our field ob- 
servations in Refs. [J, one knows that the timing and 
magnitude of phytoplankton blooms vary significantly on 
inter annual to longer timescales. For example, in some 
years, e.g., 2002, the spring bloom phytoplankton oc- 
curred in April, coinciding with vernal flowering on land, 
whereas in other years, e.g., 1999, the bloom occurred in 
February. Although there are only a handful of dynamics 
states that phytoplankton populations may exhibit stable 
equilibrium, deterministic extinction, stable population 
oscillations, and irregular fluctuations, there is a long list 
of factors that may interact to determine the community 
dynamics, e.g, competition, predation, parasitism, mutu- 
alism, age, stage and genetic structure, spatial structure 
of the habitat, climate, physical and chemical param- 
eters. Thus, a key step in analyzing the community's 
dynamics is to untangle these mixtures of interacting 
factors and to identify their essential for the observed 
dynamics. Especially, from the experimental plankton 



communities [6] [rotifer- algal] one knows that: first, the 
age structure of the predator populations is necessary to 
generate qualitatively correct predictions of population 
dynamics such as stability versus cycles. Second, the 
rapid evolution of population, including the alga and ro- 
tifers, is regarded as critical processes occurring on the 
same time scale in the ecological dynamics, and learned 
that microcosms may not just serve as a means to check 
model's assumptions, but that the results of microcosm 
studies can lead to novel insights into the function of bi- 
ological communities. From Refs. 0, BQ, we know that 
the interannual variation in zooplankton and phytoplank- 
ton species might be the result of changes in climate, and 
the photosynthetic of the phytoplankton growth strongly 
depends on the intensity of the light. Hence, these period 
factors are regarded as the external periodic forces in the 
plankton systems. The variances of species evolving in 
time and space may be well understood. 

Besides these periodic factors, there are many 
other stochastically factors causing phytoplankton- 
zooplankton blooms to transit extinction. For ex- 
ample, the effects arise from rivers' pollution on the 
phytoplankton-zooplankton ecosystem is one of the ways 
that humans affect the marine ecosystem. Large rivers 
are major mechanisms for nutrient delivery to the ocean. 
Hence, the quality of river water affects freshwater 
ecosystems and oceanic food webs. Away from this, the 
long term climatic variation also shows the stochastic fac- 
tor for the ecological system [see [10] for giving a review 
about this points]. In the phytoplankton-zooplankton 
systems, the random fluctuations also are undeniably 
arising from either environmental variability or internal 
species. To quantify the relationship between fluctua- 
tions and species' concentration with spatial degrees of 
freedom, the consideration of these fluctuations supposes 
to deal with noisy quantities whose variance might at 
times be a sizable fraction of their mean levels. For ex- 
ample, the birth and death processes of individuals are in- 
trinsically fluctuations [ll| , which become especially pro- 
nounced when the amount of individuals is small. The 
interaction between the oceanic zooplankton and fish, 
which are far from being uniformly distributed, also in- 
troduce randomness [12]. In addition, these processes 
can be regarded as parameter fluctuates irregularly with 
spaces and time. For example, phytoplankton production 
is effected by iron concentration, which can be elevated 
in surface water after rain [13]. These unavoidable fluctu- 
ations can also interact with the system's nonlinearities 
to render counterintuitive behavior, in which an increase 
in the noise level produces a more regular behavior p^ . 
Now, it is natural to ask what are spatiotemporal behav- 
iors of the spatially extended plankton system if these ex- 
ternal periodic forces and the irregular fluctuation factors 
work on. Especially, can the periodic oscillation retain to 
appear? In order to understand the stochastically force, 
we give the short introduction of the noise application in 
the physics and biology. 

Recently, research interest has shifted to the effects 
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of noise in spatially extended system [13, [lH. Well- 
known examples in zero-dimensional (ODE) systems are 
noise- induced transition [l^, stochastic resonance [13, 
Us, T9] , enhanced spatial synchronization [20] , and noise- 
sustained oscillations [2lj, [22|, [2^ [the reader can read 
an overview written by F Sagues in Ref. [14]]. More re- 
cently, effects of noise in spatially distributed systems, 
including noise-induced pattern formation [14, 24, 25], 
noise-induced fronts 0, etc, have been extended 
study by many communities. In these and other noise- 
related phenomena, multiplicative noise, which couples 
to the system state plays a very special role. However, 
prominent effect has been also found for additive noise. 
Such influence has been observed in noise-induced pat- 
tern formation [I8, 25, 2I, [H, [s^. A recent report \31\ 
demonstrated that additive noise which globally alter- 
nates between with two different monostable excitable 
dynamics yields pattern formation. Meanwhile, in re- 
cent years several theoretical investigations have been 
done on noise- induced effects in population dynamics, 
see [13, [33I, [13, di, S, [13, S, m. Noise, such as that 
from the natural variability, are inevitably presented in 
these type of systems, but its effects have not yet been 
addressed completely. In this paper, we report reso- 
nant patterns and frequency locked oscillations induced 
by additive noise and external forces in phytoplankton- 
zooplankton systems, and take into account the interplay 
among noise, external forces, and diffusion processes. In 



the following, we first give the spatially extended model 
and method we used, and then the describe our results. 



II. MODEL 

From the recent perspective of Pascual [i^ , the phys- 
ical environments play an important role in the biota. If 
the climate sharply changes, the population abundance 
will change. Especially, the spatiotemporal dynamics of 
global population abundance, such as aggregated over 
the whole space, can be approximated by mean- field- 
type equations in which the functional forms specifying 
the growth rates and interactions have been modified as 
power functions. And, the effect of interactions on local 
(or individual scales) can be represented implicitly by the 
changed form of the functions to describe interactions on 
global scales. We want to focus on the effects of a pe- 
riodic varying and the stochastic fluctuation factors to 
the spatially extended phytoplankton and zooplankton 
model. Following Scheffer's minimal approach [f4l| and 
the previous analysis [42|, |43|, |4^, 145|, we study a two- 
variable phytoplankton and zooplankton spatial model 
involving time-periodic forces and the fluctuation term 
to describe the influence of spatial noisea spatial on pat- 
tern formation. The spatially extended model is written 
as 



dp ap 

m='''^'-''^-TTVp 



h + Asin{ujt) + dpV'^p^ 



(la) 



dh CiP . in A / \ / \ 7 7 

^ h- mh - f ^ + A sm(cc;t) + r]{r, t) + dh\/ h, 



dt l^bp 



(lb) 



where the parameters r, a, b, m, n, dp, dh^ and / re- 
fer to the previous works [42|, |43| on the dimensionless 
model ([1]) absence of periodic force and noise term. Here 
p{x, t) and h{x, t) are scalar fields representing the 
concentrations of phytoplankton species and zooplank- 
ton species in the two-dimensional space. It must also 
be kept mind that the individuals of the population have 
not been at the same point in space at previous times, 
so they are a function of time. The periodic force is as- 
sumed to be sinusoidal with amplitude A and frequency 
00. The periodic force is considered as additive version 
for reason that toxins produced by different phytoplank- 
ton species has a significant role in shaping the dynam- 
ical behavior of marine plankton ecosystems [46]. The 
zooplankton population tries to avoid the areas where 
the concentration of phytoplankton is very large. The 
reason may be either dense concentration or the effect 



of toxic substance released by phytoplankton, of course 
including the others factor as introduced in previous sec- 
tion, such as human action, nature sharp changing-red 
tides and localized outbreaks and occur in coastal water 
and fronts [47]. Buskey and Stockwell [48] have shown 
in their field study that macro- and meso-zooplankton 
populations are reduced during the blooms of Chryso- 
phyte (Aureococcus anophagefferens) . These factors are 
not directly correlated with species internal parameters, 
so that the Asm{ujt) be added in the equation instead 
of ^4^(1 — p) sm{u;t). Such consideration is meaningful in 
the biology systems (see Ref. [H, EH). In Eq. ([1]), the 
stochastic factors are taken into account as the term, 
r]{r,t). Of course, due to the coupling, the noise in 
Eq. (jlbp will have its influence on Eq. (p!aj) as well. The 
noise term r]{r,t) is introduced additively in space and 
time, which is the Ornstein-Uhlenbech process that obeys 
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the fohowing stochastic partial diflFerential equation [50| : 
^ = ->,*) + >,t), (2) 

where ^{r^t) is a Gaussian white noise with zero mean 
and correlation, 

{ar,t)ar',t')) = 2sS{r-r')5{t-t'). (3) 

The colored noise rjir^t)^ which is temporally correlated 
and white in space, satisfies 

{r?(r, t)r?(r', t')) = ^ exp ( - <^(r- - r'), (4) 

where r controls the temporal correlation, and e mea- 
sures the noise intensity. 

In this work, we rely on numerical simulations of the 
model of Eqs. ([T]) and ([2]). We here consider spatiotem- 
poral evolution of this system with space white noise and 
colored noise evolving in time when the system lies within 
the regime of self-sustained Hopf oscillation. Note that 
the explicit form of the noise term represents only on 
the fiuctuating recruitment rate (or death rate) has been 
studied by Malchow et al [11], in which the parameter m 
only is regarded as noise form, i.e., m = mo + ^(t), (,{t) 
denotes Gaussian white noise. For the absence of the 
external periodic forces and colored noise, Hopf instabil- 
ity occurs and spatially homogeneous oscillation comes 
up when the parameter / less than the critical value 
fn = 0.3398 [see Appendix A] whose value depends on 
the other parameters. Moreover, from previous analy- 
sis and numerical simulations [5l|, [52[ show that spiral 
waves structure exists in the two-dimensional space and 
the spatiotemporal chaos will emerge through its far-field 
breakup when the parameter / within the mixed domain 
of the Hopf- Turing instability, but the parameters are 
outside the Turing instability domain in present paper 
[see Appendix A for a simple analysis]. In fact, from re- 
cent perspective that the 2:1 resonance phenomena on 
the spiral wave may will appear when the system un- 
dergos saddle-node or Hopf bifurcation In other 

words, the temporal period of the bifurcation patterns 
is twice the period of the primary spiral, i.e., the Hopf 
frequency ujh needs to be in a 2:1 resonance with the 
rotation frequency of the spiral wave. However, here we 
show that the system also exhibits spatiotemporal chaos 
patterns and the 4:1 and 1:1 resonances with the period of 
the external forces when fiuctuations and periodic forces 
are considered. Except when it explicitly is pointed out, 
we take parameters r = 5, a = 5, 6 = 5, m = 0.6, 
n = 0.4, / = 0.3, dp = 0.05, and dh = 0.5 throughout 
this paper. These parameters estimate by Medvinsky et 
al [42I, 'S, 'S] . From these References, we know that these 
parameters are meaningful from the ecological point of 
view. The noise intensity e and correlation time r are 
adjusted as control parameters. 



III. RESULTS 

Extensive testing was performed through numerical 
simulations to the described model ([1]) and the results 
are shown this section. In simulation, zero-fiux boundary 
conditions are used and time step At = 0.05 time unit. 
The space step Ax = Ay = 1 length unit and the grid 
sizes in the evolutional simulations are N xN {N = 200). 
The Fourier transform method is used for the determin- 
istic part in Eq. ([1]). On the discrete square lattices, the 
stochastic partial differential Eq. ([2]) is integrated nu- 
merically by applying the Euler method. Several differ- 
ent discrete methods (simple Euler, Runge-Kutta, and 
Fourier transform) was checked, and the results indicate 
that Fourier transform accurately approximates solutions 
of Eq. ([1]). On the other hand, the Fourier method offers 
the speed advantage over other numerical methods. We 
find that on the PC computer the Fourier method runs 
about 3-4 times faster than the Euler integration using 
the same time step and space step. Code is implemented 
in Matlab 7.3 and the Jft2, Jftshift, ijft2, and ifftshift func- 
tions were used for the main numerical integration. 

Although the noisy fiuctuations may sometimes caus- 
ing the variable {p and h) to be less than zero, it will 
lead to the diffusion-reaction system with cutoff effect at 
low densities when the species extinction is taken into 
account explicitly [H^]. According to the spatially ex- 
tended model ([1]), at each position in space, whenever 
the population densities fall below certain prescribed 
value e they are set to zero or sufficient small positive 
constant [Hj, [55[. From the biological point of view, 
in this paper we set that they equal to 0.0001 when 
the variables change to negative. Note that we are not 
much concerned here with the "exact" value of £, for the 
reason that an attempt to estimate the "exact" value 
would hardly make any ecological sense in terms of very 
schematic model ([1]). To compare with the numerical re- 
sults under the different cases, we used the same initial 
conditions that are randomly perturbed (the perturba- 
tions are space-independent) by homogenous equilibrium 
(p*, /i*) = (0.3944, 1.7998), except when they are explic- 
itly pointed out. 



A. Dynamics under the presence of noise only 

In the ecological systems, noise-sustained and noise- 
induced spatial pattern formations have been discussed 
in recent years [33, H, [sl, [H, [s^, [H, [H . For the conve- 
nience of discussion about the resonant pattern formation 
induced by additive noise and the external forces, here we 
first present a brief description about the affection arising 
from the fluctuations in the system ([T]). From our results 
show that the sole noise also plays a constructive role 
in the model ([1]) by maintaining (or eliminate) the large 
spatial clusters [cf. Fig. [3l [H and partly movie- l^movie- 
2^ movie-3] and enhancing the aperiodic oscillation [see 
Fig. [5]. Some typical snapshots of the spatial patterns at 



t = 300 presented in Fig. O and Fig. O [H where before 
and after the noise is added, respectively [see caption of 
Figures]. The ratio's distribution of two species is a key 
measurement in the ecological field. To characterize this 
measurement in the system ([T]), we introduce the param- 
eter, 6 (the 6 G [— 27r,27r]) and define as following: 



arg((ps,fe - P*) + i{hs,k - /^*)), 



(5) 



where 5,/c = l,2,---,A/'. ^ indicates the relation between 
the two variables of the system, called the phase angle. 
Note that it also depicts the variance of the amplitude 
relevant the two species in the space, for the reason that 
its variables are the difference around the stable points, 
and /i*. Similarly, the following parameter, £, depicts 
the average amplitude for the system in the space. How- 
ever, we find that it is equivalent to the mean value of 
the variables in this system [cf. Fig. l2l7l8l and Fig. [5], so 
we only use it in this section. 



(6) 



where s,/c = l,2,---,A^. The spatiotemporal chaotic spi- 
ral patterns appears when the noises are free [see Fig|2], 
and the large clusters can be recognized in the space. The 
jC or < p > (the < p > denotes the spatial average of p 
rather than the time average) oscillates with small fiuc- 
tuations (quasi-period) [5^], as shown in FiglSJd) (and 
Fig. E]), in which the abscissa is time and the ordinate is 
variable, C. The behavior of the system undergoes dras- 
tic changes when the noises are turned on (cf. Fig. [21 [3l 
andSj). First, the large spatiotemporal chaotic clusters 
die out gradually with the increasing of noise intensity; 
second, the oscillations of C become more obvious when 
noise intensity is within a certain regions, but it stays at 
a fixed value again if the noise intensity is further increas- 
ing. In sum, compared with Fig. [21 [3l and [H it should 
be noted that the noise plays a constructive role in the 
relation between phytoplankton and zooplankton in the 
space. 



B. Dynamics under the presence of the external 
forces only 
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FIG. 2: A typical grey-scaled snapshots of spatiotemporal 
chaotic patterns at t = 300, when the system ([T]) evolves in 
time without noise term, (a) and (b) the spatial patterns of p 
and /i, respectively; (c) the phase angle, (see the definition in 
text); (d) the amplitude, C evolves in time, [see the Movie-1, 
and additional movies format available from Ref. [531 
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FIG. 3: A typical grey-scaled snapshots of spatiotemporal 
chaotic patterns at t = 300, when the system ([T]) evolves 
in time with noise term. The meanings of the images are 
correspondent to the Fig. [2 with the values of the parameters 
are t — 1 and e — 0.001. [see the Movie- 2, and additional 
movies format available from Ref. [531 



From the previous analysis [8], it is known that the 
plankton system ([1]) has a very complex dynamical be- 
haviors such as the scenario of bifurcations and chaos 
when the system was not considered the spatial degree 
of free. However, we here consider the case that the pa- 
rameter values are within the Hopf oscillation interval 
and the spatial degree of free was taken into account. 

The parameter values under we take in previous Sec- 
tion II, and also the system is noise free in this subsection, 
thus the spatially homogeneous oscillation does not re- 
spond to the external periodic force if the amplitude A is 
below a threshold Ac, whose value depends on the exter- 
nal period Tin = 77- Above the threshold, the system ^ 



produces oscillations about period Tout with respect to 
external period Tin, this phenomena are called frequency 
locking or resonant response, i.e., when the system pro- 
duces one spike within each of the M (M = 1, 2, 3, • • • ) 
periods of the external force, that is, M : 1 resonant 
response. In present paper, the output period Tout is 
defined as follows: Ti is the time interval between the 
ith spike and the i + 1th spike, q spikes are taken into 
account and the average value of them is Tout^ where 

Tout = ^M- 

Two types of frequency locking phenomena are ob- 
served in the spatially extended model ([1]). They are 
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which plots the space-time figure for the homogeneous 
patterns to compare with the external force, and illus- 
trates that the patterns are well 4 : 1 frequency lock- 
ing. Moreover, the 4 : 1 frequency locking was de- 
mosntrated with a power spectrum as shown in Fig. [TOl 
i.e.. Tout : Tin = ^ : ^ = 4 : 1. Note that the green 
1 lines in Fig. [6] only denote the exact temporal period of 
5 the external forces rather than its values to comparing. 
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FIG. 4: A typical grey-scaled snapshots of spatiotemporal 
chaotic patterns at t = 300. Same situation as the Fig. [3] 
but s = 0.05. [see the Movie- 3, and additional movies format 
available from Ref. [57]] 
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FIG. 5: Directly compare with the mean values of variable p 
without and with noise term. Same situation as the Fig. [2]|4l 



4 : 1 and 1 : 1 frequency locking. In Fig. [61(a) and (b), 
we have plotted the temporal evolution of the variable 
p when the amplitude A is above the threshold, where 
the 4 : 1 and 1 : 1 frequency lockings take place, re- 
spectively. In these cases, the spatially homogeneous os- 
cillation patterns depend on the initial conditions, and 
two-phase patterns with a phase shift and separated by 
stationary Ising front [In particular, for the initial con- 
ditions coming from the Ref. [43] that is p(x,y,0) = p*, 
and h{x^y^O) = -\- ex -\- here e and S are parame- 
ters.], or alternative homogeneous oscillations come out 
[see Fig. El^c)]. From the random initial conditions pre- 
pared by randomly perturbing around the homogeneous 
steady state we obtain 4 : 1 resonant homo- 

geneous oscillations, as an example shown in Fig. El^c), 




600 
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FIG. 6: (color) Sequences of the mean concentration (p), the 
noise free, (a) The 4 : 1 frequency locking oscillation with 
the values of the parameters are A = 0.1 and a; = 27r x 10~^. 
(b) The 1 : 1 frequency locking oscillation with the values of 
the parameters are A = 0.1, but a; = tt x 10""*^. (c) Space- 
time diagram (right panel) displays the time evolution of the 
patterns in the one- dimensional space. For comparison, the 
curve at the left panel shows the periodic external force. In 
panel (c), the left and right window share the same ordinate 
which is time. Note that the green lines only denote the exact 
temporal period of the external forces rather than its values 
to comparing. 
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C. Dynamics under both noise and external forces 

Now, we turn on the additive noise and periodic forces 
in the system ([1]) . To check their effects on the system ([1]) 
within the 4 : 1 frequency locking regimes, we adjust 
the noise strength and correlation time. At first series 
of simulations, we adopt r = 1.0, and adjust the noise 
intensity, £. Considering the influences of noise and spa- 
tial distribution on the system, we performed simulation 
starting at e = 0.001, and then e is increased in small 
steps = 0.001 until the noise intensity is enough large 
values. The noise drastically changes the previous scenar- 
ios in the system including the spatial patterns and fre- 
quency locking when the noise intensity is strong enough 
[see Fig. [7]. First, the spatial patterns with homogeneous 
oscillation are replaced by the spatial heterogeneous os- 
cillations. This means that the distribution of the species 
may appear spatial patterns and its bloom is periodic in 
the space when some fluctuations work on the system 
[see Fig. [Tl^a)]. Second, the frequency locking (4 : 1) 
shifts into the other type frequency locking (1:1) [see 
Fig. [lb)]. 

From our numerical results [cf. Fig. [7], one could con- 
clude that homogeneous oscillations depends on the noise 
intensity. The oscillations are only slightly perturbed, or 
not affected when the noise intensity e is small. As e is in- 
creased, the homogeneous oscillating patters lose its sta- 
bility and associate with the frequency locking changed. 
Figure [TJa) (middle column) shows that oscillating two- 
phase pattern with local spotted patterns will appear at 
£ = 0.001. When the noise intensity further is increased, 
the spotted patterns show aggregation behaviour, and 
become a large clusters in the space. Figure [^a) (right- 
hand column) shows a clear spatial oscillation patterns 
with a large clusters will appear at e = 0.050. In Fig.fTJb) 
we have plotted the temporal evolution of the variable (p) 
when the noise intensity is taken different values. It is 
interesting to realize that the resonant patterns and the 
frequency locking do not alter by increasing the noise in- 
tensity within the large interval [e.g. 0.0005 to 0.001; 
0.05 to 0.10]. The variable {h) also exhibits the similar 
characteristics in model. 

In order to elucidate the evolutional processes of the 
spatial structure, we have depicted typical spatiotempo- 
ral patterns of the system ^ within one period Tout^ as 
shown in Fig. [8] and the movie- 6. The other periods also 
exhibit the similar characteristics. 

The previous Figs. [6] and [3 exhibit an aspect about the 
frequency locking and resonant pattern which have re- 
ceived considerable attention in the recent years, namely, 
the response of the system to a periodic forces may be 
enhanced by the presence of noise [13, E, [l^, HO, lE^, 113] • 

The role of temporal correlation r of the colored noise 
is the other significant except it intensity in inducing and 
controlling the spatial patterns formation and transition 
of the resonant patterns. Now, it is natural to ask what 
is the effect consequence of the temporal correlation of 
the colored noise. Especially, the phase diagram of the 
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FIG. 7: Spatial frequency locking patterns and the time se- 
ries of the mean value (p) of the concentration p{x,y) with 
respect to different noise intensities s. (a) Grey-scaled snap- 
shots of spatial frequency locking patterns with different noise 
intensities s, where the noise intensity is 0.001 and 0.05 for 
middle column and right-hand column, (b) The change of the 
frequency locking with respect to the different noise intensi- 
ties. Same situation as the Fig. [6{a) but e. [see the Movie-4, 
Movie-5 for Fig.[7ja), and additional movies format available 
from Ref [52]] 

e — r parameter space. In order to well understand the 
phase transition by the influence of temporal correlation 
r, we performed a series of simulations, fixing the r and 
scanning the noise intensity, £, when the frequency lock- 
ing evidently changes in the long term [here we run the 
time up to 5000] and recorded the data. Figure [9] sum- 
marizes the results from the numerical simulations, in 
which the region A and B are corresponding to the 1:1 
frequency locking and 4:1 frequency locking respectively. 
Figure [9] depicts the transition point of frequency locking 
shifts toward higher values of the noise intensity as the 
correlation time is increased, that is, r softens the effect 
of the noise. 



IV. CONCLUSION AND DISCUSSION 

It is known that, an external periodic forces applied to 
a nonlinear pendulum can cause the pendulum to be- 
come entrained at a frequency which is rationally re- 
lated to the applied frequency, a phenomenon known as 
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FIG. 8: Typical spatial pattern formation to the forces-noise 
system ([!]) within the 4:1 frequency locking region as in case 
Fig. [Z{b) [e = 0.005]. (a) Grey-scaled snapshots of spatial 
patterns of the variable p and the phase angle, 0. [the time 
from the left to right], (b) shows the time series of the mean 
concentration p corresponding to the snapshots in (a) . [see the 
Movie-6 for Fig.[8ja), and additional movies format available 
from Ref. [57]. Note that the movie is 25 frames per second]. 
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FIG. 9: Phase diagram in e — r parameter space, with A = 0.1 
and a; = 27r X 10~^. There are two states: 4:1 region (B) and 
1:1 region (A). The solid lines is a least-square fit of circles 
data. 



the frequency- locking [6l|. A recent theoretical analy- 
sis shows that an array of coupled nonlinear oscillators 
can exhibit spatial reorganization when they subjected to 
external periodic forces [62j]. In this paper, we present a 
spatial version of phytoplankton-zooplankton model that 
includes some important factors such as external periodic 
forces, random fluctuations, and diffusion processes. The 
modified system based on the original model by Schef- 
fer [41, 42], but here it can exhibit frequency- locking phe- 
nomena and resonant patterns. Our results show that 
the noise and the external periodic forces play a con- 
structive role in the Scheffer's model: first, the noise can 
enhance the oscillation of the species' density and format 
a large clusters in the space. Second, the external peri- 



odic forces can induce 4:1 and 1:1 frequency- locking and 
spatially homogeneous oscillation to appear. Third, the 
resonant patterns are observed in the system when the 
spatial noise and external periodic forces are turned on, 
moreover the frequency-lockings transit when the noise 
intensity increased. These result from outside the do- 
main of Turing instability. In the Appendix [Aj we fur- 
ther provide a simple Turing linear stability analysis with 
the help of the numerical calculation by using the Maple 
software. 

It is worth emphasizing that the frequency lockings and 
resonant patterns only appear when noise and the exter- 
nal forces present in the model ([1]). We here consider 
the cases when the unforced system lies outside the Tur- 
ing instability [see Appendix [A] for details] . The typical 
power spectra P{oo^) related to the density of the phyto- 
plankton for the oscillations within 4:1 frequency locking 
is analysis, as shown in Fig. [TOl The power spectra re- 
lated to the density of zooplankton [not shown here] is 
similar. The presence of a prominent and well defined 
peak in the power spectrum of Fig. [10] at a nonzero fre- 
quency characterize oscillating behavior. From Fig. [TOl 
one could see that the oscillation behavior relates to the 
external periodic forces and the natural frequency, coh 
[see the Appendix [X] of the system. Significantly, oscil- 
lations are enhanced in the system when the single noise 
term presents. These results indicate that the oceanic 
plankton bloom may partly due to the interplay between 
external forces and stochastic factors instead of determin- 
istic factors. Our results also may help us to understand 
the effects arising from undeniable subject to random 
fluctuations in oceanic plankton bloom. 

In this paper, we study the effects of noise and external 
forces on spatially extended phytoplankton-zooplankton 
system in static media (no prominent the advection 
term). However, in the oceanic ecological systems the 
biological processes among the species are in a fluid 
environment, such as turbulent flows or chaotic advec- 
tion. Recently, a few authors consider the mixing of the 
flow and diffusion processes using a well-known standard 
model of chaotic advection [3, 23, 63, 64, 65, £6, 67,^] in 
the excitable media, especially, Reigada et al [g^ found 
that plankton blooms also can be induced by turbulent 
flows. So, a further study in our work will analyze the 
effects arising from the mixing of the advection, diffu- 
sion, the external forces and noise. It is also interesting 
to apply the problem of oceanic plankton bloom and the 
spatial structure observed by the field as a introduction 
in Section I. 
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APPENDIX A: STABILITY ANALYSIS WITH 
THE HELP OF MAPLE 



Now we consider the stability of this positive equilib- 
rium. The Jacobian matrix at the positive equilibrium 
{p\h*) is 



J 



dp 



dp 



ip*,h*) 



dgi 
dh 

dg2 



dh 



ip^,h*) 



(A4) 



J(p*,/i*) 



TABLE I: The eigenvalues of the Jacobian matrix ()A4I 



/ 0.30 


0.3397 


0.3398 


0.3399 


Ai,2 .0473±.822H 


.0001±.8040z 


.0000±.8039z 


-.0001=b.8039z 



Its eigenvalues, Ai,2, were listed in the Table [J with the 
different parameter values /. The first type of instability 
is associated with a Hopf bifurcation in the spatially uni- 
form system. This instability will appear at the thresh- 
old fn and oscillate homogeneously with the natural fre- 
quency loh = Ini(Ai^2)- From Table H one could see that 
the natural frequency equals to 0.8221 if / = 0.30. 

Now we consider the spatial inhomogeneous perturba- 
tion by using the techniques of Koch and Meinhardt . 

Sp = Spo ex.p{X't + ikr), Sh = Sho exp(A t + ikr)^ (A5) 

with wave number vector k = (kx^ky), \^Po\, \^ho\ <C 
1 and imaginary unit i'^ = —1. Due to the zero-fiux 
boundary conditions, kx and ky take only discrete values 



k^ = /c^ = n7r/L,n = 0,1,2 •• 



(A6) 



This appendix is devoted to the numerical analysis of 
unforced system (jAip . 



dp X «P . 

— - = rpil — p) —a, 

dt l^bp ' 



dh CLP . in ^h'^ 
-a — ma — j- 



dt l^bp 



(Ala) 



(Alb) 



Setting the left-hand sides of system ()Aip to zero, we 
obtain 



92{p, h) 



gi{p, h) = rp{l -p)- -r^K (A2) 
1 + op 

ap 



l^hp 



h — mh — f 



ap 
bp' 



(A3) 



The Eqs. (jA2|) and (jA3|) exist a unique interior equilib- 
rium (p*, /i*), when parameter / less than 0.445. 



where L is the size of the space. Each k^ is associated 
with a "frequency" cj^, which can be a complex number. 
The functions conik^) are found by substituting expres- 
sion ()A5p into following equation 



dp 



-P) 



ap 



l^bp 



h + dpV'^p, 



(A7a) 



^ h-mh- f ^ , + dhV^h. (A7b) 



dt l^bp 



Retaining terms up to first order in 5p and 5h^ we get 
linearized equation: 



J' 



I ( 

6h 



0, 



(A8) 



with 
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J' = 



dp 



dg2 



dg2 
dh 



dgi 
dh 



ip*,h*) 



ip*,h*) 

-dhik^y 



(A9) 



TABLE II: The eigenvalues of the Jacobian matrix ()A9|) 



/ 


0.30 


0.3124 


0.3125 


^1,2 


-0.0704zb0.8177z 


-0.0003±0.8174z 


0.0003=b0.8174z 




0.6539 


0.3480 


0.3444 



The perturbation ampHtudes Spo and 5ho can be different 
from zero if and only if the det J' = 0. 

Turing instability is expected to occur for finite /c^ > 0, 
and at least one of the real parts becoming greater than 
zero. If / is taken as a control parameter, the critical 
point is reached if the determinant J' has a pure imagi- 



nary eigenvalue, i.e., / = /t, Ke{X[ 2) = 0- And for the 
critical wave number can be written as 

zdpdh 

Jii and J22 are the elements of the Matrix J. 

Table [IT] lists the eigenvalues of the Jacobian ma- 
trix ()A9p and the critical wave number /c^, from which 
one can see that our simulation outside the Turing insta- 
bility regime when the / = 0.30. 

The Maple program available on request. 
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